Chaotic traveling rolls in Rayleigh-Benard convection 
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^ : Abstract 

In this paper we investigate two-dimensional (2D) Rayleigh-Benard convection using direct nu- 
merical simulation in Boussinesq fluids of Prandtl number P = 6.8 confined between thermally 
conducting plates. We show through the simulation that in a small range of reduced Rayleigh 
| number r (770 < r < 890) the 2D rolls move chaotically in a direction normal to the roll axis. 

^f) ■ The lateral shift of the rolls may lead to global flow reversal of the convective motion. The chaotic 



traveling rolls are observed in simulation with free-slip as well as no-slip boundary conditions on 



o 

the velocity field. We show that the traveling rolls and the flow reversal are due to an interplay 



between the real and imaginary parts of the critical modes. 
03 1 PACS numbers: 47.20.Bp, 47.20.-k 



The Rayleigh-Benard convection (RBC) a a a is an extensively studied system for 
investigating a range of interesting phenomena like instabilities jj], pattern-formation jjj, 
chaos jg, ill, spatio-temporal chaos, and turbulence Q, 0, [l0|. The convective flow is charac- 
terized by two non-dimensional numbers: the Rayleigh number R and the Prandtl number 
P. Two-dimensional (2D) stationary rolls 2| appear as primary instability in Boussinesq 
fluids confined between thermally conducting boundaries when the Rayleigh number R is 
raised just above a critical value R c . The convective dynamics at higher values of R depends 
on Prandl number. Interesting convective dynamics are observed when the Rayleigh number 

nnnnn □ 

is increased beyond R c |4|, [ll|, [12j, ll3j, [14[ . Busse [12| found traveling waves along the axis 
of the cylindrical rolls as secondary instability in low Prandtl number fluids that makes the 
convection three-dimensional (3D). Traveling waves are also found in rotating RBC in cylin- 



drical geometry 151 ] . Another kind of traveling waves is known to occur in two-dimensional 
(2D) RBC in binary mixtures 16, l^j], where the straight rolls move in a direction perpen- 



dicular to the roll axis. Such behaviour is not reported in two-dimensional RBC in pure 
fluids, although many interesting features have been investigated (e.g., 0,0,0). 

In this paper we present results of two-dimensional numerical simulation of the RBC 
in a pure Boussinesq fluid of Prandtl number P = 6.8. We have used free-slip as well as 
no-slip boundary conditions on the velocity field at the horizontal plates. In the horizontal 
direction, we have assumed periodic boundary conditions on all the fields. We report our 
results as a function of reduced Rayleigh number r = R/R c . For free-slip boundaries, at 
r fa 130 we observe a transition from one oscillatory state with only real or imaginary 
values of critical modes W101 and #101 to another oscillatory state where both the real and 
imaginary parts are non-zero. At later r, near 770 < r < 890 we find chaotic traveling waves 
in a direction normal to the roll-axis. The sudden lateral shift in the roll system leads to 
either jitters in the convective flow or flow reversal of the convective motion. The chaotic 
travelling waves are also observed in the direct numerical simulation (DNS) of the RBC with 
no-slip boundaries. We find links between the chaotic traveling waves or rolls and the flow 
reversal. 

We consider 2D convection in an extended layer of Boussinesq fluid with thermal ex- 
pansion coefficient a, kinematic viscosity u, thermal diffusivity k that is enclosed between 
two flat conducting plates separated by distance d and heated from below. The adverse 
temperature gradient is AT/d, where AT is the temperature difference imposed across the 




FIG. 1: (a) Timeseries of K(Wioi), 3r(Wioi), 3?(^ioi) 5 9(#ioi) at reduced Rayleigh number r = 
700 and Prandtl number P = 6.8. The timeseries has two leading frequencies, and it shows 
quasiperiodic behaviour, (b) The projection of the phase space in K(Wioi) — ^(^101) that elucidates 
the quasiperiodic behaviour of the system. 

bounding plates. The nondimensional equations are 

3w 

— + (u ■ V)u = -Va + RP9z + PV 2 u, 

— + (u-v)o = u ?J + v 2 e. 

at 

Here the Rayleigh number R = Qg ^p d is the ratio of the buoyancy and the dissipative 
forces, while the Prandtl number P = u/n is the ratio of thermal diffusive time 7#, = d 2 /n 
and the viscous diffusive time r vis = d 2 /v. The temperature perturbation 9 due to the 
convective flow vanishes at thermally conducting boundaries. The realistic no-slip boundary 
conditions imply that the velocity field v = (u,0,w) = at the boundaries. The idealized 
free-slip boundaries imply d z u = w = at the flat plates. The value of R c is 277r 4 /4 pa 657.5 
for free-slip boundaries, and R c ~ 1707.8 for no-slip boundaries. 

We have carried out our simulations of the RBC in 2D using pseudo-spectral method [Lj 
with stress-free boundaries. The grid resolution used for the simulations with the free-slip 
boundaries is 256 x 256 with an aspect ratio of 2v^2 : I. The reduced Rayleigh number 
r = R/R c is varied from 1.01 to 10 3 . We have employed fourth-order Runge-Kutta (RK4) 
scheme for time stepping with the time step varying from 1 x 10~ 4 to 1 x 10 -6 in thermal 
diffusive time units. The code was validated using the results reported by Thual [rsl ]. 

To solve no-slip convective flow, we employed finite-difference two-step precedure with 
the Adam-Bashforth-Crank-Nicolson scheme. The first step involves pressure calculation 
using Poisson's equation, and the second step deals with time-advancing using fourth-order 



central explicit scheme with enhanced spectral resolution. We computed the no-slip 
flow for r = 830 and P = 6.8 and compared the results with that obtained with free-slip 
simulation. 

First we report the results of 2D simulation with free-slip boundary conditions. We 
observe long time behaviour of the critical 2D modes Wioi and #101. Here the three subscript 
denote the wavenumber indices of the Fourier modes along x, y, and z directions, where z 
is the vertical direction. In our simulation the second subscript is zero. We start our 
simulation with real Fourier modes. For reduced Rayleigh number r < 80, the system shows 
time independent steady convection. For 80 < r < 130, we observe simple oscillations in 
all the modes. Here all the Fourier modes remain real. Similar observations have been 
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made previously by others Us Il9l . |20j. As r is raised further, the critical Fourier modes 
(Wioi and 9 W i) become complex. For 130 < r < 145, the imaginary parts of W W i and 9 W i 
become nonzero and they oscillate around a zero mean, while the real parts oscillate around 
a nonzero mean. 

For 145 < r < 770 the imaginary parts of the complex modes also oscillate around a 
finite mean. This is an oscillatory instability, which makes the critical modes complex. 
The real and imaginary parts of a Fourier mode interact with each other through nonlinear 
interaction with higher order modes. When the simulation is started with purely imaginary 
modes (except #oon, which is always real due to the reality condition), then the roles of real 
and imaginary parts of all complex modes get interchanged. 

The oscillations of both real and imaginary parts become strongly anharmonic for 660 < 
r < 770. This leads to temporally quasiperiodic flow. Fig. [T](a) illustrates the real and 
the imaginary parts of the critical mode Wiqi and d±oi- These modes show two leading 
frequencies. The projection of phase space on Wiqi — #101 shows quasiperiodic flow as seen 
in Fig. [IJb). Similar behavior has been observed in experiments on RBC by Swinney and 

Gdiub y. 

The convection becomes chaotic for 770 < r < 890. Figure [2] displays the real and the 
imaginary parts of the critical modes W^ioi and 8 W i for r = 830. The projection of the phase 
space on $l(Wioi) — ^t(0 W i) plane is illustrated in Fig. [3]that indicates the chaotic nature of 
the system. 

In the chaotic regime, the real and imaginary parts of the critical modes interact strongly 
through nonlinear coupling with higher order Fourier modes. As illustrated in Fig. [2], the 



250 r 



1 



-250 



la) 




time (d Ik) 



FIG. 2: Timeseries of the real and imaginary parts of the critical velocity mode W±ox ( a )j an d 
the temperature mode #101 (b) for r = 830. The origin on dimensionless time axis is suitably 
chosen after the steady-state has been reached. The timeseries of the modes shows jitters (small 
fluctuations) and flips. The fluctuations in the temperature mode is stronger than in that in the 
velocity mode. 
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FIG. 3: Projection of the phase space in the 3i(Wioi) — 3i(0ioi) plane at r = 830. The figure 
illustrates the chaotic behaviour of the system. 



Fourier modes show jitters (small fluctuations) as well as change in sign. The velocity 
mode W\ox an d temperature mode #101 appear to be in approximate phase with each other. 
However the fluctuations in the temperature mode near the jitter or the flip is stronger than 
the corresponding fluctuations in the velocity mode. The interval between two flips is of the 
order of one thermal diffusion time scale, while the duration of the flips is much shorter. 

The physical interpretation of the chaotic timeseries (Fig. [2]) yields interesting insights 
into the traveling wave instability and flow reversal. We provide a simple argument to explain 
this. If the critical modes were the only modes in the system, then the vertical velocity in 
real space would be 3?(Wioi exp(ik c x) sin(rwr2:)). When Wioi changes from a real value A to a 




FIG. 4: Travelling rolls in two-dimensional Rayleigh-Benard convection at r = 830. The top panel 
of the figures is for the simulation with free-slip boundaries and the bottom panel is for the no-slip 
boundaries. The temperature field is represented in color with blue being the coldest region. The 
arrows represent the velocity field. The three images at the top panel show the field configurations 
at three time. At t = 1.2514 the central region has hot fluid that shifts leftward subsequently. By 
t = 1.35 the whole roll pattern is shifted approximately by half a wavelength. Similar features are 
seen for no-slip convective flow, except that the flow pattern moves right. 

complex value Aexp(i(f)), it would lead to a shift of the roll pattern in horizontal distance by 
4>/k c . This is the observed traveling wave instability near the bifurcation. Since the temporal 
change in the phase of the critical mode is chaotic, the movement of the rolls would also 
be irregular. During the flip of the critical Fourier modes, ~ n, which corresponds to a 
horizontal shift of A/2. Under this situation, the vertical velocity reverses its direction at a 
given location that leads to global flow reversal in the system. The simulated flow contains 
many modes apart from the critical modes, yet the large-scale ones play a dominant role. 
Hence, the bifurcation of the large-scale Fourier modes from real to complex play a major 
role in initiating traveling wave instability and flow reversal. 

We illustrate the above dynamics using Fig. H] that contains two panels of images: the top 
ones are for the free-slip boundaries, while the bottom ones are for the no-slip boundaries. 
The two panels of Fig. H] illustrate three snap-shots of the velocity and temperature at three 
different instances of time. The top panel (free-slip) illustrates the flow fields at t = 1.2514, 



t = 1.28, and t = 1.35 respectively for a different run. The temperature fields are shown 
in color coding with blue being the coldest region, while the velocity fields are shown using 
arrows. At t — 1.2514, the hot fluid is rising near the center. The roll pattern moves 
leftward as time advances. At time t = 1.35, a cold fluid parcel is falling down near the 
center. Thus the top panel illustrates the traveling roll in RB convection under free-slip 
boundary condition. In the bottom panel, similar feature is observed for no-slip boundary 
condition however for a shorter time. Here the rolls have moved only a smaller fraction 
of a wavelength. Note that the rolls are traveling perpendicular to the roll axis along the 
periodic direction in both the results, and the interval between two consecutive flow reversal 
is of the order of one diffusive time unit (refer to Fig. [2]). 

Figure H] also depicts flow reversal phenomenon. If we place a probe near the center of 
the system, we would observe the fluid to be hot and moving upward at initial time. At the 
end, the above mentioned probe would measure colder fluid flowing downward. Thus a flow 
reversal due to traveling rolls is clearly illustrated. 



We have created two movies to illustrate the traveling rolls and flow reversa 
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by taking 



23] . The 



many frames between the reversal. These movies can be seen in our website 
first movie 22J depicts the traveling rolls near a flow reversal regime. By the end of the 



movie (in approximately 0.1 thermal diffusive time unit), the roll has shifted approximately 
by half a wavelength. Similar features are observed in the second movie 23J] that depicts 
the convective flow with no-slip boundaries; here the roll pattern moves to the right. 

We have performed several three-dimensional (3D) simulations for Prandtl number 6.8, 
and a range of Rayleigh number in the chaotic and turbulent regime. We observe the above 
features of traveling rolls and flow reversal in 3D as well; these reults will be reported in due 
course. 
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arge-scale flows has been observed in many experiments and numerical simu- 
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281 ] . Some of the experiments have been performed for large Rayleigh 



numbers (of the order of 10 9 to 10 12 ). The physics of the above phenomenon has been of 
great interest in the recent past. Our results are for two dimensional RBC at moderately 
low Rayleigh numbers. Yet, the simple mechanism of traveling rolls in the periodic direction 
and its consequence to the phenomenon of flow reversal may be of interest in more complex 
situations. The large Rayleigh number regimes are under investigation. 

To summarize, we observe succession of patterns in two-dimensional RBC under free- 



slip boundaries. Our focus has been on the chaotic traveling rolls observed in the range of 
770 < r < 890. We show that the generation of traveling rolls is due to the bifurcation 
of critical modes from real values to complex values. This kind of traveling roll is possible 
along the periodic direction that may be realized in practice along the azimuthal direction 
in a cylindrical container. It is interesting to note that several traveling roll states and flow 
reversals have been observed in cylindrical geometry. Also, the chaotic traveling rolls lead 
to flow reversal when the amplitude of the critical Fourier modes switches sign. 
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